AD-A244  474 

liililiiii 


.  .hiu 

-(g) 


IMPROVED  PROPAGATION  MODELS 
IRREGULAR  MEDIA 


FINAL  REPORT 


by  Charles  L.  Rino 


30  June,  1991 


U.  S.  ARMY  RESEARCH  OFFICE 


Grant  Number:  DAAL03-89-C-0011 


Vista  Research,  Inc. 

100  View  Street  P.  0.  Box  998 
Mountain  View,  CA  94042 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


VISTA  RESEARCH,  INC. 

100  View  Street  •  P.  D.  Box  SS8 
Mountain  View,  CA  34042  ■  [41 5]  SSB-II^I 


92-00756 


05 


149 


IMPROVED  PROPAGATION  MODELS 
IRREGULAR  MEDIA 

FINAL  REPORT 

by  Charles  L.  Rino 

30  June,  1991 

U.  S.  ARMY  RESEARCH  OFFICE 

Grant  Number:  DAAL03-89-C-0011 

Vista  Research,  Inc. 

100  View  Street  P.  O.  Box  998 
Mountain  View,  CA  94042 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THE  VIEW,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS  REPORT 
ARE  THOSE  OF  THE  AUTHOR(S)  AND  SHOULD  NOT  BE  CONSTRUED  AS  AN 
OFFICIAL  DEPARTMENT  OF  THE  ARMY  POSITION,  POLICY,  OR  DECISION, 
UNLESS  SO  DESIGNATED  BY  OTHER  DOCUMENTATION. 


REPORT  DOCUMENVAIIUN  PAGE 


OMB  No.  0?04-0iea 


Puoiic  'WJOrttnq  ouraert  *Of  thu  cctiection  oi  informacion  eMimxted  to  i  no«r  o^t  f^%oot'%€.  meiuotn^  th«  time  tof  reviewing  instructions,  seercnmg  exsting  data  sources, 

gathering  and  maintaining  the  data  needed,  and  ccmoieting  and  reviewing  the  collection  ol  information,  ^end  comments  regarding  tnis  burden  estimate  or  any  other  asoect  of  tn»s 
collection  ol  mtormanon.  including  suggestions  tor  reducing  this  ouroen  to  Washington  Headquarters  Services.  Directorate  tor  information  Ooeractons  arsd  heoorts.  12!5  ietlerson 
Dans  Highway.  Suite  1204  Arlington,  vA  22202-4)02.  and  to  me  Office  of  Management  and  Oudgct.  PsperworK  fteduaion  Proiea(07C4>0l88).  Washington.  DC  20503 


1.  ACENCY  USE  ONLY  (Leive  bUnk) 


2.  REPORT  DATE 

30  June,  1991 


4.  TITLE  AND  SUBTITLE 

"IMPROVED  PROPAGATION  MODELS  IRREGULAR  MEDIA' 


6.  AUTHOR(S) 


S.  FUNDING  NUMBERS 
U.S.  Army  Research  Office 
Grant  Number: 

DAAL03-89-C-e«  00|| 


by  Charles  L.  Rino 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS<ES) 

Vista  Research,  Inc. 

P.O.  Box  998 

100  View  St.,  Ste.  #202 

Mountain  View,  CA  94042 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AOORESS{ES) 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


Project  1024 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


Altt)  <3.4^  yd.  A-G-S'-S 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Miximum  200  words) 


ABSTRACT 


This  final  report  suinmarizes  the  most  important  results  obtained  from  the  re¬ 
search  supported  under  the  Army  Research  Office  Grant  No.  DAAL03-89-C-0011. 
We  have  completed  work  that  was  initiated  under  a  previous  ARO  contract  (DAAL03- 
87-COOl)  to  develop  a  consistent  and  rigorously  correct  formalism  that  can  be  applied 
to  both  discrete  and  continuous  random  media.  That  work  has  led  to  a  new  ana¬ 
lytic/computational  method,  which  we  call  the  mutual  interaction  method  (MIM).  The 
method  has  been  successfully  applied  to  problems  involving  a  discrete  scatterer  near  a 
random  surface.  It  also  forms  the  basis  for  highly  efficient  computational  algorithms 
for  numerical  simulations. 


14.  SUBJEa  TERMS 

Propogation 
Random  media 


IS.  NUMBER  OF  PAGES 

13  (including  cover) 


1«.  PRICE  COOE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


MSN  7S40-01 -280-5500 


18.  SECURITY  CLASSIFICATION 


UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF  ABSTRACT 


8 


Standard  Form  298  (Rev  2-39) 

^rncriMd  by  ANSI  Std  239*’8 
298-102 


GENERAL  IMSTRUCTIONS  FOR  COMPLETIWG  SF 


The  Report  Documentation  Page  (RDP)  is  used  in  announcing  and  cataloging  reports.  It  is  important 
that  this  information  be  consistent  with  the  rest  of  the  report,  particularly  the  cover  and  title  page. 
Instructions  for  filling  in  each  block  of  the  form  follow.  It  is  important  to  stay  within  the  lines  to  meet 
optical  scanning  requirements. 


Block  1.  Agency  Use  Only  fLeave  b/ank). 

Block  2.  Report  Date.  Full  publication  date 
including  day,  month,  and  year,  if  available  (e.g.  1 
Jan  88).  Must  cite  at  least  the  year. 


Block  12a.  Distribution/Availabilitv  Statement. 
Denotes  public  availability  or  limitations.  Cite  any 
availability  to  the  public.  Enter  additional 
limitations  or  special  markings  in  all  capitals  (e.g. 
NOFORN,  REL,  ITAR). 


Blocks.  Type  of  Report  and  Dates  Covered. 
State  whether  report  is  interim,  final,  etc.  If 
applicable,  enter  inclusive  report  dates  (e.g.  10 
Jun87-30Jun88). 


Block  4.  Title  and  Subtitle.  A  title  is  taken  from 
the  part  of  the  report  that  provides  the  most 
meaningful  and  complete  information.  When  a 
report  is  prepared  in  more  than  one  volume, 
repeat  the  primary  title,  add  volume  number,  and 
include  subtitle  for  the  specific  volume.  On 
classified  documents  enter  the  title  classification 
in  parentheses. 


Blocks.  Funding  Numbers.  To  include  contract 
and  grant  numbers;  may  include  program 
element  number(s),  project  number(s),  task 
number(s),  and  work  unit  number(s).  Use  the 
following  labels: 


C  -  Contract 
G  -  Grant 
PE  -  Program 
Element 


PR  -  Project 
TA  -  Task 
WU  -  Work  Unit 

Accession  No. 


DOD  '  See  DoDD  5230.24, 'Distribution 
Statements  on  Technical 
Documents.' 

DOE  '  See  authorities. 

NASA  -  See  Handbook  NHB  2200.2. 

NTIS  '  Leave  blank. 


Block  12b.  Distribution  Code. 

DOD  •  Leave  blank. 

DOE  '  Enter  DOE  distribution  categories 
from  the  Standard  Distribution  for 
Unclassified  Scientific  and  Technical 
Reports. 

NASA  •  Leave  blank. 

NTIS  •  Leave  blank. 

Block  13.  Abstract.  Include  a  brief  (Max/mum 
200  words)  factual  summary  of  the  most 
significant  information  contained  in  the  report. 


Block  6.  Author(s).  Name(s)  of  person(s) 
responsible  for  writing  the  report,  performing 
the  research,  or  credited  with  the  content  of  the 
report.  If  editor  or  compiler,  this  should  follow 
the  name(s). 

Block?.  Performing  Organization  Name(s)  and 
Address(es).  Self-explanatory. 

Block  8.  Performing  Organization  Report 
Number.  Enter  the  unique  alphanumeric  report 
number(s)  assigned  by  the  organization 
performing  the  report. 

Block  9.  Soonsoring/Monitorino  Agency  Name(s) 
and  Address(es).  Self-explanatory. 

Block  10.  Sponsorino/Monitoring  Agency 
Report  Number.  (If  known) 

Block  11.  Supplementary  Notes.  Enter 
information  not  included  elsewhere  such  as; 
Prepared  in  cooperation  with...;  Trans,  of...;  To  be 
published  in....  When  a  report  is  revised,  include 
a  statement  whether  the  new  report  supersedes 
or  supplements  the  older  report. 


Block  14.  Subject  Terms.  Keywords  or  phrases 
identifying  major  subjects  in  the  report. 

Block  15.  Number  of  Pages.  Enter  the  total 
number  of  pages. 

Block  16.  Price  Code.  Enter  appropriate  price 
code  (NTIS  only). 

Blocks  17.  *19.  Security  Classifications.  Self- 
explanatory.  Enter  U.S.  Security  Classification  in 
accordance  with  U.S.  Security  Regulations  (i.e., 
UNCLASSIFIED).  If  form  contains  classified 
information,  stamp  classification  on  the  top  and 
bottom  of  the  page. 

Block  20.  Limitation  of  Abstract.  This  Nock  must 
be  completed  to  assign  a  limitation  to  the 
abstract.  Enter  either  UL  (unlimited)  or  SAR  (same 
as  report).  An  entry  in  this  block  is  necessary  if 
the  abstract  is  to  be  limited.  If  blank,  the  abstract 
is  assumed  to  be  unlimited. 


9 


Standard  Form  298  Back  (Rev.  2-89) 


I  statement  of  the  Problem 


In  continuous  random  media,  one  invariably  uses  the  parabolic  approximation  to  the 
wave  equation.  Thus,  the  development  of  moment  equations  that  characterize  the 
random  field  proceeds  from  a  model  that  excludes  a  priori  wide-angle  scattering  and 
backscatter.  While  attempts  have  been  made  to  rectify  both  limitations,  the  formu¬ 
lations  used  are  intractable  or  inconsistent.  It  is  desirable  to  use  a  formulation  that 
accommodates  backscatter  and  wideangle  scatter  at  the  outset. 

In  discrete  random  media,  the  formalism  developed  by  Flody,  Lax,  and  Twersky 
is  more  often  used.  The  problem  development  is  setup  so  that  a  self-consistent  com¬ 
putation  of  the  complete  multiple  scattering  interactions  amoung  all  the  scatters  is 
accommodated.  It  is  known  that  self-consistent  interaction  computations  can  be  set  up 
as  solutions  to  differential  equations  as  well  as  implicit  summations  of  all  interactions 
(diagram  methods).  Whereas  the  continuous  media  problem  generally  proceed  from  a 
system  of  restricted  differential  equations,  the  discrete  problem  more  often  proceeds  an 
exact  diagram  system.  It  is  desirable  to  use  a  common  formulation  that  preservs  the 
self-consistent  interaction  fields  but  can  be  transformed  to  diagram  form. 

We  believe  that  our  current  work  has  achieved  this  goal  and  provides  some  new 
insights  into  the  structure  of  scattering  problems. 

II  Summary  of  Principal  Results 

Completing  work  that  was  initiated  under  a  previous  ARO  contract  (DAAL03-87-C- 
002),  we  have  developed  a  consistent  formalism  akin  to  the  multiple-phtise-screen  model 
that  is  rigorously  correct  for  both  discrete  and  continuous  random  media.  For  example, 
the  mutual  interaction  form  of  the  propagation  equations  for  continuous  random  media 
can  be  written  as  follows: 


aV'+CKiz) 

dz 

ik  ft  ^  -  dK' 

0) 

d^~{K\z) 

dz 

(2) 

The  integral  term  invol’ 

ves  the  total  wave  field 

xi>{K\z)=  ^+(K;z)-|-  ^~[K\z), 

(3) 

and  8t{K\z)  is  the  spatial  Fourier  transform  of  the  relative  permittivity  variation. 
Insofar  as  we  know,  these  equations  have  not  been  used  previously. 
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In  [1],  we  describe  the  results  of  a  detailed  study  of  the  conditions  under  which  a 
hirearchy  of  moment  equations  can  be  derived  from  (1)  and  (2).  We  show  that  two 
distinct  approximations  approximations  are  involved.  The  first  approximation  requires 
small  local  perturbations  as  a  sufficient  condition  for  evaluating  the  functional  deriva¬ 
tives  that  appear  in  the  Novikov-Furutsu  (NF)  theorem.  This  assumption  alone  permits 
the  development  of  a  closed  hirearchy  of  differential  equations  for  the  signal  moments. 
The  second  approximation  requires  that  over  short  distances  the  spectral  components 
propagate  like  plane  waves  in  a  homogeneous  medium.  Under  the  first  and  second  ap¬ 
proximations,  the  most  general  form  of  the  integral  equations  for  the  second  order-signal 
moments  that  depend  only  on  the  spectral  density  function  of  the  permittivity  fluctu¬ 
ations  can  be  developed.  Upon  simplifying  these  equations  by  using  the  narrow-angle 
scatter  approximation,  we  reconfirmed  an  earlier  finding,  namely  that  backscatter  en¬ 
hancements  are  negligibly  small  in  continuous  random  media  under  the  usual  conditions 
of  the  Markov  approximation. 

In  earlier  work  [2],  we  performed  a  comparative  analysis  of  solutions  to  the  one¬ 
dimensional  form  of  the  wave  equation.  The  method  of  invariant  imbedding  provides  a 
rigorous  alternative  solution  at  the  level  of  the  first  approximation.  Thus,  the  second 
approximation  can  be  tested.  We  found  that  for  the  first-order  moments  of  the  field, 
the  Markov  solution  shows  little  error;  however,  the  total  flux  within  the  medium  ad¬ 
mits  a  systematically  increasing  error  as  the  portion  of  the  incident  flux  that  has  been 
backscattered  incre<ises.  The  invariant  imbedding  method  does  not  resolve  the  internal 
flux  into  directed  wave  fields,  but  for  a  lossless  slab  of  sufficient  depth,  all  of  the  inci¬ 
dent  intensity  will  be  backscattered.  The  breakdown  of  the  second  approximation  is  of 
practical  concern  for  sound  and  EM  propagation  in  highly  inhomogeneous  media  such 
as  soil  or  mud.  Similar  limitations  most  likely  will  apply  in  nonsparse  discrete  random 
media. 

Unfortunately,  configurational  averaging  in  discrete  random  media  has  not  produced 
rigorous  solutions  at  the  level  of  the  invariant  imbedding  solution.  Thus,  in  [3],  we 
assumed  that  the  fields  at  the  boundary  planes  of  the  nth  slab  are  uncorrelated  with 
the  scatterers  within  the  slab.  With  this  assumption,  it  is  possible  to  derive  a  closed 
hierarchy  moment  equations  for  the  vector  fields.  We  showed  that  the  characteristic 
equation  for  the  mean  field  propagating  in  a  statistically  homogeneous  medium  has 
four  eigenvalues  corresponding  to  orthogonally  polarized  waves  propagating  in  the  ±2 
directions.  For  the  scalar  wave  equation,  the  paired  eigenvalues  are  given  as 

A  =  ±  (a  ±  ^(1  -h  a++){\  +a—)-  a+-a-+  -h  A^  )  ,  (4) 

where 

A  =  (a++-a--)/2,  (5) 
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and  represents  the  ensemble  average  of  the  corresponding  scattering  function  for 
single  slab  but  normalized  by  the  slab  thickness.  When  =  a  =  and  o-'*'”  = 
O’"*"  =  a**,  the  result  agrees  with  a  well-known  result  derived  by  Twersky,  which  is 
commonly  used  to  estimate  the  extinction  effects  of  backscatter. 

The  uncorrelated  field  hypothesis  appears  to  be  valid  for  sparse  media,  but  its 
limitations  presently  are  uncertain.  Thus,  we  have  not  yet  extended  the  computations 
to  higher-order  moments. 

We  have  applied  the  MIM  method  to  practical  problems  involving  an  object  scat¬ 
tering  near  a  highly  rough  surface.  Under  a  separate  contract,  we  applied  the  method 
to  investigate  the  hypothesis  that  enhanced  acoustic  surface  reverberation  in  the  ocean 
is  caused  by  subsurface  bubble  clouds  [4].  These  computations  provide  a  good  example 
of  the  comparative  ease  with  which  the  results  from  numerical  simulations  and  known 
scattering  functions  can  be  combined  via  MIM.  A  surface-scatter  simulation  was  de¬ 
veloped  for  dynamically  evolving  nonlinear  ocean  surface  realizations  as  one  input  to 
the  MIM  computation  [5].  The  bubble  cloud  was  modeled  as  a  cylindrical  void  with 
known  scattering  characteristics,  although  analytic  continuation  must  be  used  to  extend 
the  Bessel  series  to  accommodate  evanescent  waves.  The  general  problem  of  a  particle 
scattering  near  a  rough  surface  is  reviewed  in  [6]. 

The  scattering  problem  is  simplified  considerably  when  the  integration  implied  by 
the  @  operator  can  be  replaced  by  a  simple  product.  In  [7],  we  describe  a  class  of  scat¬ 
tering  objects  for  which  the  scattering  operator  can  be  evaluated  algebraically.  This 
is  important,  because  the  general  scattering  interaction  of  a  wave  field  with  a  particle 
requires  the  equivalent  of  an  integral  operation  whether  formulated  in  the  spatial  or 
Fourier  domains.  The  T-matrix  method  effectively  converts  the  integral  to  an  infinite 
series  of  spherical-harmonic  modes.  A  finite  subset  of  the  coefficients  are  manipulated 
with  the  T-matrix  approach.  For  objects  that  do  not  deviate  strongly  from  spheros- 
ity,  this  method  is  efficient;  however,  for  a  broader  class  of  scattering  objects,  direct 
computation  of  the  scattering  function  as  described  in  [8]  may  be  more  efficient. 

Finally,  although  it  was  not  formally  part  of  our  current  ARO  contract,  we  used  our 
surface-scatter  code  to  study  the  statistics  of  backscatter  enhancements.  Our  results 
confirmed  earlier  findings  that  speckle  statistics  remain  gaussian  distributed  through 
the  region  of  the  backscatter  enhancement  [9].  We  have  also  simulated  backscatter 
enhancements  for  highly  nongaussian  surfaces  at  near-grazing  incidence.  Results  from 
our  simulation  codes  have  been  made  available  to  other  ARO  researchers. 


Some  unpublished  work  as  yet  unpublished  is  described  in  more  detail  in  the  next 
section. 
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Ill  The  Relation  of  MIM  to  Order  N  Algorithms 

III.l  Background 

Wayne  Chew  [10]  has  described  a  recursive  algorithm  that  evidently  can  achieve 
computational  efficiency  for  a  system  of  N  mutually  interacting  scatterers.  The  scheme 
has  been  compared  to  0{N^)  method-of-moments  solvers  to  verify  that  it  is  indeed 
solving  comparable  problems  with  substantially  improved  efficiency  [llj.  More  recently, 
an  even  faster  implementation  of  the  algorithm  has  been  described  [12].  The  latter 
scheme  is  structurally  similar  to  the  mutual  interaction  method  [3],  which  we  believe 
to  be  capable  of  similar  efficiencies.  The  relation  is  difficult  to  see,  however,  because 
Chew’s  method  uses  the  T-matrix  formalism;  however.  Chew’s  method  can  be  developed 
more  generally. 

Consider  a  system  of  N  isotropic  scatterers.  The  structure  of  the  solution  will  show 
how  the  computational  efficiency  is  achieved  and  how  to  reformulate  the  method  in 
terms  of  other  scattering  operators  [7].  For  an  incident  wave  iP{t)  scattering  from  a 
system  of  N  scatterers,  the  total  field  at  any  exterior  point  r  can  be  written  as 

V'^(r)  =  ^(r)  +  H  (6) 

i=i 

wliere  i/i?  represents  the  field  scattered  by  the  jth  scatterer  in  the  presence  of  all  other 
scatterers.  This  field  can  be  related  linearly  to  the  total  wave  field  incident  upon  the 
scatterer,  as  follows: 

V' ;(r)  =  JJJ^^  5j(r,  r')  dr\  (7) 

where  5j(r,  r')  depends  only  on  the  attributes  of  the  jth  scatterer  and  the  integration 
is  carried  out  over  its  volume  Vj.  The  total  field  incident  upon  the  jth  scatterer  admits 
the  representation 

=  .^(r)  +  E  ///,  S»(r,  r')«(p')  dr'.  (8) 

n^j 


The  system  of  N  equations  implied  by  (8)  can,  in  principle,  be  solved  for  the  un¬ 
known  fields  which  need  be  determined  only  within  the  confines  of  each  scatterer 

where  5'j(r,  r')  is  finite.  The  formal  solution  to  (8)  can  be  written  as 

=  Tiif,  (9) 


S 


where  is  the  jth  row  of  the  N  x  N  matrix  operator  and  is  a  column  vector 
whose  nth  element  consists  of  the  incident  field  evaluated  in  the  volume  of  the  nth 
scatterer.  In  terms  of  lyv,  (6)  can  be  written  as 


(10) 


As  a  practical  matter,  the  problem  is  solved  when  the  fields  are  determined.  As 
with  any  linear  system  of  equations  more  computation  is  required  to  determine  lyy 
itself;  moreover,  the  scattered  field  must  be  evaluated  in  the  region  of  interest.  This 
formalism  is  attributed  to  Foldy,  Lax,  and  Twersky  (FLT). 

III. 2  Isotropic  Scatterers 

An  object  located  at  is  called  a  point-like  isotropic  scatterer  if  its  scattering  function 
Sj  has  the  following  property: 

///  Si(r,r>j,(r')<ir'=  (11) 

In  (11),  the  right  hand  side  is  the  scattered  wave  observed  at  r,  0  is  an  arbitrary 
incident  wave,  hj  is  a  constant  characterizing  the  object’s  scattering  strength,  and  G  is 
the  outgoing  Green’s  function  defined  by 

_  _  exp{ifclr-r'|} 

=  '  4^|r-?r 

Consider  a  rightward  traveling  incident  wave  V>(r)  multiply  scattered  by  N  point-like 
isotropic  scatterers.  The  FLT  equations  simplify  to 

^^(r)  =  V’(r)  +  Y.  Tj)  (13) 

i=i 


and 


—  V’(rj)  +  Y  (14) 

n«l 

n/j 
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wlieie  Gjn  =  G(rj,r„).  The  latter  equation  can  be  written  in  matrix  form  as 


1 

—  h2G\2 

~^3f7i3  .  . 

-hf^GiN 

V*! 

—h^G2\ 

1 

-h^G-a  ■ 

—  hf^'G2N 

'f'N 

V’2 

—hiGsi 

“^2^32 

1  .  . 

—hffGsN 

= 

V’a 

-hiG^i 

—  h2Gs2 

.  1 

1 

Jn. 

.  . 

If  we  let  Xw  represent  the  inverse  of  the  matrix  in  (15),  the  formal  solution  is 

(16) 

We  could  equally  well  work  in  terms  of  the  scattered  fields  Either  way, 

the  direct  solution  to  this  system  oi  equations  would  require  0{N^)  operations 


III.3  A  Fast  Algorithm  for  Isotropic  Scatterers 


Now  consider  a  system  of  N  isotropic  scatterers  with  individual  strength  hj  as  described 
in  Section  III.  We  will  formulate  the  solution  in  terms  of  which  represents  the  total 
scattered  field  at  the  location  of  the  jth  particle,  Tj.  In  terms  of  4'^,  the  total  scattered 
field  can  be  written  as 


Recall  that 

(16) 

where  is  the  total  field  incident  upon  the  jth  particle.  Moreover,  V’(r)  is  a  freely 
propagating  wave  field.  Thus,  if  the  field  is  known  on  a  plane,  say  a,t  z  =  zq,  it  can 
be  transformed  to  any  other  plane,  say  at  z  =  Zj,  by  using  the  propagation  operator 
denoted  as  P:0. 

This  property  can  be  used  to  develop  an  efficient  recursive  algorithm  for  computing 
For  a  single  isolated  particle,  If  two  particles  are  present,  it  follows 

that 


V'2  =  ht  [V»i  +  Gi2V’|]  (19) 

=  h2  [^2  +  (20) 
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The  terms  in  square  brackets  represent  the  direct  and  scattered  wave  impinging  upon 
the  the  particular  scatterer.  For  the  propagating  wave  field, 


rf>j  =  PjorpQ. 


(21) 


The  direct  parallel  to  Chew’s  [10]  development  can  be  seen  by  the  following  equivalences: 


fj 

=  rki 

(22) 

^i(i) 

=  h, 

(23) 

=  G., 

(24) 

PjO 

=  P.0 

(25) 

a 

=  V’o 

(26) 

We  now  substitute  (19)  into  (20)  and  use  the  relation  Pi 0^*0= A 2 ^20^*0  to  obtain 


/ 2  I  [1  +  P l2G2lh\]  „ 


(27) 


If  we  define 

T  —  A2^*21^i] 

~  1  -  /li/t2Gl2G2i  ’ 


(28) 


we  see  that  the  scattered  field  at  the  second  scatterer  is  given  in  terms  of  a  constant 
multiplying  the  incident  field  at  the  location  of  the  scatterer,  p2oV’o-  Substituting  (28) 
into  (19)  and  manipulating  the  incident  field  yields  the  complementary  relation  for  the 
first  scatterer,  namely 

P|(2)  =  ^l[l  +  G\2T2(2)P 2i]-  (29) 


Because  Tj(\)  =  /i,,  we  see  that  the  2-particle  T  elements  are  given  in  terms  of  the  single 
particle  T  elements. 

For  an  n-particle  system. 


V’n  =  =  Pj(„)PjoV’0- 


From  (17),  it  follows  that 


V7"+' 


= -/■ + ct; + E 


+!• 


J=1 


The  form  of  (31)  together  with  (30)  implies  that 

V’n+l  —  [PjoV’O  +  f^j.n+lV’n+l]  > 


(30) 


(31) 


(32) 


8 


which  is  the  general  form  of  (19).  Also,  for  an  n  +  1-particle  system, 


V’n+l  —  ^n+l  l^^n+1,0^0  +  ^<Jn+l,iV’n+l  >  (33) 

which  is  the  generalized  form  of  (20).  Proceeding  exactly  as  with  the  two-particle 
system,  we  substitute  (32)  into  (33)  and  manipulate  the  result  to  obtain 


=  h  ^  ~1~  IZj=i  ^n+\jTj(n) Pj, n+l  „  . 

Vn+I  ^n+1  1  1  'T  ^n+l.oV’Oj 

*  Xnjljj.n+l 


from  which  it  follows  that 


^  d-  (^n+i,jTj^n)Pj,n+l 


-*n+l(n+l)  -‘n+l(l).  _  ,  C  T'  C  ’  V'*®/ 

1  “n+1  Z^j=l  ^n+\,j-i  j(n)^j,n+l 

Substituting  (35)  into  (20)  and  manipulating  the  incident  field  completes  the  recursion: 

^i(n+l)  =  7j(n)  [l  +  G^i,n+l^n+l(n+l)^n+l,j]  •  (36) 


Equations  (35)  and  (36)  arc  functionally  identical  to  Chew’s  [10]  equations  (17)  and 
(18).  Although  Pjo  is  an  operator  that  would  generally  require  an  FFT  to  evaluate, 
the  total  number  of  operations  required  to  evaluate  (35)  and  (36)  is  0{N^).  From  this 
discussion  it  is  clear  that  the  computational  efficiency  arises  by  virture  of  the  fact  that 
the  incident  wavefield  is  freely  propagating  and  therefore  can  be  mapped  from  position 
to  position.  Tne  scattered  fields  have  the  same  property  outside  the  scatterers.  The 
T-matrix  addition  theorem  propagates  the  incident  fields  (via  /3ij)  and  the  scattered 
fields  (via  a,j)  from  the  location  of  one  particle  to  another. 

We  have  used  plane  and  spherical- wave  propagators,  but  the  procedure  is  general 
More  important,  there  are  other  ways  to  exploit  and  improve  on  these  efficiencies  within 
the  context  of  the  general  operator  formalism  described  in  [8].  As  Chew  has  noted, 
the  summation  terms  in  (35)  and  (36)  can  be  put  into  a  common  form,  which  is  the 
aggregate  scattering  function  for  the  collection  of  particles.  In  terms  of  the  aggregate 
scallering  function,  the  solution  is  effectively  a  succession  of  two-element  scattering 
systems.  We  have  noted  that  the  Mutial-Interaction-Equations  could  be  solved  in  this 
fashion,  with  improved  efficiency.  Chew’s  results  show  that  this  conjecture  is  correct. 


III. 4  An  Example 

For  the  simple  case  of  an  incident  plane  wave  (a  delta-function  in  the  spatial  Fourier 
domain),  the  propagation  operator  is  multiplicative,  whereby  all  operations  are  purely 
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algebraic.  In  the  T-matrix  formalism,  this  is  an  M  =  1  system.  We  have  programmed 
the  recursive  algorithm  and  compared  the  results  to  the  exact  solution  for  both  accuracy 
and  efficiency.  With  the  caveat  that  the  system  is  not  close  to  being  singular,  in  which 
case  the  denominator  in  (35)  and  (36)  is  zero,  the  errors  are  negligible.  The  CPU 
time  for  a  SUN  SPARC  2  computer  is  shown  in  the  figure.  The  upper  curve  is  the 
LU  decomposition;  the  lower  curve  is  the  recursive  algorithm.  The  variations  from  the 
expected  and  curves  is  most  likely  due  to  the  dynamic  operating  environment. 
The  main  point  is  that  for  the  large  system,  the  computation  time  is  reduce  from  more 
than  one  hour  to  a  few  minutes. 


Figure  1:  CPU  time  comparison  of  direct  and  recursive  solutions. 
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